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Abstract 

We calculate Coulomb breakup of the neutron rich nucleus on a Pb 
target at the beam energy of 244 MeV/nucleon within the framework of a 
hnite range distorted wave Born approximation theory that is extended to 
include the effects of projectile deformation. In this theory, the breakup 
amplitude involves the full wave function of the projectile ground state. 
Calculations have been carried out for the total one-neutron removal cross 
section (cr_i„), the neutron-core relative energy spectrum, the parallel mo¬ 
mentum distribution of the core fragment, the valence neutron angular, and 
energy-angular distributions. The calculated cr_i„ has been compared with 
the recently measured data to put constraints on the spin parity, and the 
one-neutron separation energy (S'n) of the ^'^Mg ground state (^’^Mg^s). The 
dependence of cr_in on the deformation of this state has also been investi¬ 
gated. While a spin parity assignment of 7/2“ for the is ruled out 

by our study, neither of the 3/2“ and 1/2’'' assignments can be clearly ex¬ 
cluded. Using the spectroscopic factor of one for both the 3/2“ and 1/2’'' 
conhgurations and ignoring the projectile deformation effects, the Sn values 
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of 0.35±0.06 MeV and 0.50±0.07 MeV, respectively, are extracted for the two 
configurations. However, the extracted S'„ is strongly dependent on the spec¬ 
troscopic factor and the deformation effects of the respective conhguration. 
The narrow parallel momentum distribution of the core fragment and the 
strong forward peaking of the valence neutron angular distribution suggest a 
one-neutron halo conhguration in either of the 2 p 3/2 and 2 si /2 conhgurations 
of the ^^Mg ground state. 

Keywords: Coulomb breakup of ^^Mg, deformation effects, one-neutron 
removal cross section, relative energy spectra, parallel momentum 
distributions, angular distributions 
PACS: 24.10.Eq, 25.60.Gc, 27.30.+t 


1. Introduction 


With the advances made in the technology of producing nuclear species 
with relatively large neutron [N) to proton (Z) number ratios, it is now 
possible to extensively study nuclei near the neutron-drip line with Z > 8. 
During the last three decades measurements performed on mass, radius and 
spectroscopy of such nuclei have shown that they have structures that are 
at variance with those of their “near the line of stability” counterparts fsee. 


variance witn tnose oi tneir near tne line oi staDiiity 

•g., BBBSSSHSSHQQ, and HQHE 


171, list Il9|,l20|,l2]|). 


With the advent of new generation of radioactive ion beam facilities, it has 
now become possible not only to produce medium mass neutron rich nuclei 
in the vicinity of the magic numbers but also employ them as projectiles to 


initiate reactions (e.g., breakup) on nuclear targets [22|, l23|. This provides 


excellent opportunity to perform quantitative study of the single-particle 
structure and the shell evolution in this region. 

The notion of “magic” numbers is one of the most fundamental concepts 

,[3- If large gaps occur between groups of 


in nuclear structure physics 


single-particle orbits that are completely hlled with nucleons (neutrons or 
protons), then these nucleon numbers are called “magic”. The seven most 
established magic numbers are 2, 8, 20, 28, 50, 82, and 126. However, in 
several nuclei near the neutron-drip line, modihcations to this shell struc¬ 


ture have been observed [15|. In these cases the magic numbers evolve as a 


function of the neutron number - old magic numbers may disappear while 
new ones emerge and conventional shell gaps may break down. The region 
where abrupt onset of changes in the magic numbers appears, is called is- 
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land of inversion [2^. For example, rapid changes in nuclear structure and 


vanishing of the N = 8, and 20 shell gaps have been seen in neutron rich 
nuclei ^^Be [^, and ^^Mg Q and respectively. Examples of 


21 and 4"Si 28 


N = 28 shell quenching have been observed in 
It is suggested in Ref. [2^ that island of inversion near N = 20, 28 comes 
about because of the fact that the h'{sd)~^{fpY intruder conhgurations (here 
u represents a relative neutron state), in which two neutrons from the sd 
shell are excited to the fp shell, become so low in energy that they form the 
ground states for Z = 10 — 12 and iV = 20 — 22 nuclei. This suggestion was 
conhrmed subsequently by mass measurements of the neutron rich isotopes 
of Ne, Na and Mg nuclei j^. Recently, this behavior has been shown to 
be a general phenomena that should occur for most standard shell closures 
far from the line of stability, and the mechanism behind this is found to be 
related to the importance of the nucleon-nucleon tensor interaction 3^. It 


is obvious that due to the intruder states, the single particle structure of the 
ground states of nuclei lying within island of inversion will not be the same 
as that emerging from the usual hlling of the shell model states. 

The mixing of neutron n-particle-n-hole {np — nh) intruder conhgurations 
of h'{sd)~^{fp)"‘ character to the ground state, causes large deformation to 
nuclei in island of inversion near N = 20, 28, which is conhrmed by the mea¬ 
sured low excitation e nerg ies and B{E2) values of the hrst excited states (see, 
eg.. Refs. [^, @, @, 13, 16, [m 1)- It has been emphasized j^, that 

the deformation may also account for the enhanced binding energies mani¬ 
fested in some of the known nuclei in this region. The collectiye properties 
of neutron rich nuclei near = 20 region are rather well described by state- 
of-the-art Monte-Carlo shell model calculations that allow for unrestricted 
mixing of neutron particle-hole conhgurations across the shell gap [^, 36 


In Ref. [^, nuclei in the neighborhood of neutron-drip line haye been 
systematically investigated in a model where one-particle motion is described 
within spherical as well as deformed potentials. It has been concluded in this 
work that nuclei in the region of = 20 — 28 are most likely to be deformed. 

The root mean square radius (RMS) of a deformed nucleus becomes ef¬ 
fectively larger than that of a spherical one. This enhances the total reaction 
cross section (ctr) that depends on RMS radii of the projectile and the target 
nuclei. Large interaction cross sections (which are almost the same as the 


an) have been measured for ^^Ne [39l. |40 


and ^"‘“^^Mg nuclei. This 
has led to the conclusion that the isotopes have strong deforma¬ 
tion [i^, 43, 441. For the Mg case, these studies suggest that while 
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and are spherical, 25,29,33-38]y[g more likely to be deformed. 

The discovery of halo structure in some of the drip line nuclei is another 
important progress made in the studies of nuclei with large N io Z ratio near 
the limits of nuclear stability ji^, 113,0 0,0,0, Is^. As the neutron-drip 
line is approached nuclei experience weakening of the neutron binding energy, 
which leads to some special effects. The sudden rise of interaction cross 
sections with increasing N in some of these nuclei can be attributed to the 
extended density distribution (s) of the valence neutron(s). This decoupling 
of the valence neutron(s) from the tightly bound core and the extension of 
the corresponding wave function to much larger radii have been referred to 
as neutron halo. This phenomena has been seen earlier in lighter nuclei like 
iiLi, iiRe, 0. 

In recent years, there has been considerable interest in hnding out if 
halo conhgurations also exist in nucl ei ly ing near the neutron-drip line in 


the vicinity of island of inversion ^ 1^. ^^Ne with Z = 10 and = 21, 


is a promising candidate to have a one-neutron halo configuration because 
the one-neutron separation energ y of this nucleus is quite small (0.29 ± 1.64 
MeV j^, or 0.06 ± 0.41 MeV ]^). Indeed, such a structure has been 


suggested for this isotope by Coulomb breakup studies [22|. This has been 
further supported by measurements of the interaction cross sections for 
Ne isotopes incident on a target at the beam energy of 240 MeV/nucleon, 
where it was found that for ^^Ne the interaction cross section was much larger 
than that of any other Ne isotope. Recently, measurements of for 
isotopes on target at the beam energy of 240 MeV/nucleon have been 


reported in Ref. j4l| . From a similar reasoning, it was suggested in this study 
that ^'^Mg [Z = 12, V = 25) that lies in V = 20 — 28 island of inversion, is 
also a candidate for having a one-neutron halo structure. This was reinforced 
by measurements of Coulomb breakup of ^^Mg on C and Pb targets at the 
beam energy of 244 MeV/nucleon in Ref. j^. 

The observation of the halo phenomena in the heavier nuclei lying in is¬ 
land of inversion, signals major changes in the shell evolution in these nuclei 
as compared to that seen in the spherical ones. One of the conditions for 
the halo formation is that the loosely bound neutron in the nucleus occu¬ 
pies a low orbital angular momentum state (£ = 0 or 1) in order to reduce 


the centrifugal barrier effects that prevent it from spreading out [53|. In 
fact, in well established cases of light one-neutron halo nuclei like ^^Be and 
^®C, the grou nd states have predominant s-wave neutron plus core conhgura¬ 
tions 1^, l58|, l59|, l6y . According to the conventional shell model evolution. 
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one expects to see the domination of the I/ 7/2 orbit in nuclei in the vicinity 
of = 20 — 28. This would not favor the halo formation because a larger 
centrifugal barrier would prevent the I = 3 neutrons to extend too far out in 
the space. Therefore, a signihcant contribution from the s- or p-wave orbits 
has to be there in the ground state structure of these nuclei to minimize the 
centrifugal barrier. Thus, the existence of halo structure would imply a sig¬ 
nihcant modihcation of the shell structure that involves considerable mixing 
of the intruder states like 2 p 3/2 or 2 si /2 into the ground states of these nuclei, 
which also leads to the appreciable deformation of these states. Therefore, 
the halo formation in heavier nuclei in island of inversion region has strong 
correlation with the shell evolution and the presence of deformation. 

Coulomb breakup reaction, in which the valence neutron is removed from 
the fast projectile in Coulomb helds of heavy target nuclei, provides a conve¬ 
nient tool to inve stig ate the halo structure in the neutron-drip line nuclei (see, 
e.g., the review |62|). It places constraints on their electric dipole response 

the post form hnite 
re- 


571 IGIL l62| . A class of theories of this reaction (e.g. 


range distorted wave Born approximation (FRDWBA) theory |^, [63 


quires realistic wave functions to describe the relative motion between the 
valence neutron and the core in the ground state of the projectile. Thus 
by comparing the calculations of the cross sections with the measured data, 
one can directly probe the ground state structure of the projectile. Narrow 
widths of the parallel momentum distributions of the core fragments provide 
a robust signature of the presence of halo structure in the projectile nuclei 
as they imply a larger spatial spread of the fragments in their ground states. 

Recently j^, the FRDWBA theory of Coulomb breakup reactions has 
been extended to include the deformation of the projectile by using a de¬ 
formed Woods-Saxon potential to describe the valence neutron-core relative 
motion. Hence, this provides a microscopic theoretical tool to study the 
Coulomb breakup of neutron-drip line nuclei lying in island of inversion in 
the vicinity of = 20 — 28 and to investigate the correlation between halo 
formation and the shell evolution and deformation. In the hrst application of 
this theory, the Coulomb breakup of ^^Ne on a Pb target at the beam energy 
of 234 MeV/nucleon was investigated . Comparison of calculated and the 


measured (T_i„ (of Refs. [22|, l6^) suggested that the ground state of ^^Ne 
is most likely to have a ^°Ne ® 2p3/2i^ conhguration. The value of the one- 
neutron separation energy Sn is found to be correlated to the quadrupole 
deformation parameter {(32). For (32 between 0.0 — 0.5, Sn varies between 
0.24 — 0.58 MeV. The calculated full width at half maximum (FWHM) of 
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the parallel momentum distribution of ^°Ne fragment is closer to that of the 
core fragment seen in the breakup reaction of the established halo nucleus 
^^Be. This strongly suggests that the ground state of ^^Ne has a one-neutron 
halo structure in the 2 ^ 3/2 state. 

^^Mg is the most neutron-rich bound odd-mass Mg isotope. However, the 
experimental information about its mass and the ground state spin parity is 
not available, even though mass systematics suggest that it is a very weakly 
bound system with S'„ in the range of 0.16 ± 0.68 MeV j^. Therefore, it 
is another promising candidate for having a one-neutron halo structure in 
island of inversion near = 20 — 28. However, a I/ 7/2 conhguration for 
its ground state that would result in the conventional spherical shell model, 
will suppress the halo formation due to the high centrifugal barrier. Hence, 
a signihcant modihcation of its spherical shell structure with introduction of 
the intruder conhgurations having s- and p-wave states, is necessary for this 
nucleus to have a halo like structure. This also implies that its ground state 
should be deformed. 

The aim of this paper is to investigate the one-neutron removal cross 
section of ^^Mg on a Pb target at the beam energy of 244 MeV/nucleon 
within the FRDWBA theory of Coulomb breakup reactions. By comparing 
the a-in calculated within this theory with the corresponding experimental 
data we would like to extract most plausible spin parity for the ground state 
of this nucleus. We attempt to put constraints on the large uncertainty in 
its Sn value. We investigate the effect of the ground state deformation of 
this nucleus on the values of (T_i„ deduced in our analysis. Furthermore, we 
make predictions for the observables such as relative energy spectra of the 
valence neutron-core fragments, parallel momentum distribution of the core 
fragment, and the angular distribution of the valence neutron as a function 
of deformation. Our study is expected to provide more understanding about 
the evolution of the shell structure in island of inversion from = 20 to 28, 
about which some conflicting results have been reported recently [l^ 21 


Furthermore, our study is expected to quantify the presence of a neutron 
halo structure in ^^Mg and provide information about its correlation to the 
ground state deformation of this nucleus. 

The study of Coulomb breakup of ^^Mg is also of interest in astrophysics 
because it provides an indirect way to determine the rate of the radiative 
neutron capture reaction ^®Mg(?7,, 7 )^'^Mg, which is of importance in the study 
of the r-process nucleosynthesis in supernovae (b^ . 

In the next section, we present our formalism where we recall some impor- 
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tant aspects of the FRDWBA theory of breakup reactions and its extension 
to include the deformation of the projectile ground state. The results of our 
calculations are presented in section 3, where we discuss the one-neutron re¬ 
moval cross section, the relative energy spectra of the fragments, the parallel 
momentum distribution of the core fragment, the angular and energy distri¬ 
bution of valence neutron as a function of the projectile deformation. The 
summary and conclusions of our study are presented in section 4. 

2. Formalism 

The breakup reaction of a projectile a into the core fragment b and the 
valence neutron n, in the Coulomb held of a target t can be represented as 
a + t ^ b + n + t. We assume that target nucleus remains in its ground state 
during the breakup process. Thus this is also known as the elastic breakup 
reaction. The chosen coordinate system is shown in Fig. 1. The position 



t 


Figure 1: The three-body coordinate system. The charged core fragment, the valence 
neutron and the target nucleus are denoted by b, n and t, respectively. 


vectors satisfy the following relations: 


r = Fj — ari, a = -, 

-b rrih 

fn = 7 ri + <^ri, 6 = —^—, 7 = (1 - ad). 

rub -b mt 


( 1 ) 

( 2 ) 


The starting point of the FRDWBA theory of Coulomb breakup is the 
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post-form T-matrix of the reaction given by 


T 


d^dridvixi ^*(kb, 


rn)<h;(^n)Hn(ri)^i+H^a, Tl, f*). 


( 3 ) 


The functions y are the distorted waves for the relative motions of b and n 
with respect to t and the center of mass (c.m.) of the b+t system, respectively. 
The functions $ are the internal state wave functions of the concerned parti¬ 
cles dependent on the internal coordinates The function ri, rj) is 

the exact three-body scattering wave function of the projectile with a wave 
vector ka satisfying outgoing boundary conditions. The vectors kf, and k„ 
are the Jacobi wave vectors of b and n, respectively, in the hnal channel of 
the reaction. The function Vf,„(ri) represents the interaction between b and 
n. As we concentrate only on Coulomb breakup, the function x[~^(kf,,r) is 
taken as the Coulomb distorted wave (for a point Coulomb interaction be¬ 
tween the charged core b and the target) satisfying incoming wave boundary 
conditions, and the function describing the relative motion of 

the neutron with respect to the target, is just a plane wave. It may be noted 
that within this approach the fragment-target interactions are treated to all 
orders. 

In the distorted wave Born approximation (DWBA), we write 

ri, r*) = <ha(^a, ri)x^+)(k„, r^), (4) 


The assumption inherent in Eq. (4) is that the breakup channels are very 
weakly coupled and hence this coupling needs to be treated only in the hrst 
order. We expect this approximation to be valid for those cases where there 
are no resonances in the n + b continuum. Most of the neutron halo sys¬ 
tems come in this category. For those cases where higher order effects of 
fragment-fragment interaction are non-negligible, the applicability of this 
method would be limited. Ideally, a rigorous description of the breakup 
process of all types of projectiles would require the use of Faddeev type of 
three-body methods that include Coulomb potentials in the fragment-target 
and fragment-fragment (if required) interactions. A few such calculations 
have become available recently alth oug h they are conhned to the breakup 
reactions on a proton target j^, |^, 70. Til . 

In Eq. (4) the dependence of 4>a on ri describes the relative motion of 
the fragments b and n in the ground state of the projectile. The function 
xi’''^(ka, Fj) is the Coulomb distorted scattering wave describing the relative 





motion of the c.m. of the projectile with respect to the target, satisfying 
outgoing wave boundary conditions. 

The integration over the internal coordinates ^ in the T-matrix gives 

imjfi 

The wave function can be expressed in terms of its radial and angular 

parts assuming a particular partition, in which the relative motion between 
n and b has an orbital angular momentum i as 

$„(ri) = ^/rau,(rl)yf(rl), (6) 

where C‘^S is the spectroscopic factor for the given partition. In Eq. (5), i 
is coupled to the spin of n and the resultant channel spin j is coupled to the 
spin jh of the core b to yield the spin of a (jq). The partition represented by 
Eq. (6) will be retained even if the potential Vbn(ri) is deformed. 

The T-matrix can now be written as 

T = {imjnlln\jjj) {jbfj-bjfj-\jaf^a)i^il3em{^b, K] kg), (7) 

imjfi 


where 


^ ik/i.rri 


i(k;„ k„; ka) = J dridvixi ^*(k6,r)e' 

X Vbn{ri)ue{ri)YJ^{fi)xi^\K, r*)- 


( 8 ) 


with /dirn being the reduced T-matrix and with i = Y 1. 

Eq. (8) involves a six dimensional integral which makes the computa¬ 
tion of quite complicated. The problem gets further acute because the 
integrand involves a product of three scattering waves that exhibit oscilla¬ 
tory behavior asymptotically. In the past calculations have been simplihed 
by using app roximate methods, such as the zero-range ap pro ximation (see 


e.g., ( 72 I, ITSI, IT^) or the Baur-Trautmann approximation [75|, |76| that led 


to the factorization of the T-matrix into two terms; each involving three- 
dimensional integrals (we refer to for a detailed discussion). However, 
these methods are not valid for breakup reactions at higher beam energies 
and for heavier projectiles that can have non-s-wave ground states. 
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In our FRDWBA theory we use a method that was proposed in Ref. 



xt\h,r) = (9) 


Eq. (9) represents an exact Taylor series expansion abont r* if K = —iV^ is 
treated exactly. However, instead of doing this we employ a local momentnm 
approximation (LMA) where the magnitnde of momentnm K is taken to be 


K(n) = \/^(e-v(k)). 


( 10 ) 


Here m is the rednced mass of the b — t system, E is the energy of particle b 
relative to the target in the c.m. system and V{TV) is the Conlomb potential 
between b and the target separated by TZ. Thns, the magnitnde of the mo¬ 
mentnm of K is evalnated at some separation TZ, which is held hxed for all 
the valnes of r. The valne of TZ was taken to be eqnal to 10 fm. For reactions 
nnder investigation in this paper, the magnitude of K remains constant for 
distances larger than 10 fm. Due to the peripheral nature of the breakup 
reaction, the region 7^ > 10 fm contributes maximum to the cross section. 


Furthermore, the results of the calculations for these reactions, at the beam 
energies under investigation, are almost independent of the choice of the di¬ 
rection of momentum K Therefore, we have taken the directions of K 


and kb to be the same in all the calculations presented in this paper. 

Substituting Eq. (9) into Eq. (8), the reduced amplitude is obtained in a 
factorized form as 


= j dvic D)xi'^nqa, d) 

X J drie~'^-''^Vbn{ri)ui{ri)Y[^{ri), (11) 

where, Q = — aK. The hrst integral in Eq. (11), is the dynamics 

part in the Coulomb breakup and is expressed analytically in terms of the 
Bremsstrahlung integral j^. The second integral in Eq. (11) contains the 
projectile structure information. 
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We now introduce deformation in potential Vbn{Yi) in Eq. (11). Follow¬ 
ing j^, we write the axially symmetric quadrupole-deformed Woods-Saxon 
potential (without taking the spin-orbit term) as 

Vbn{^i) = 14s(n) - /32fc(ri)y2°(fi)- (12) 


We take the Woods-Saxon form for the potential 14)s(’"i), and write, V^s^ri) = 
K«s/(^i)) where is the depth of the potential and /(ri) describes its 
shape. f{ri) and k{ri) are dehned as 


fin) = 


1 -I- ’ 


kin) = RVi 


0 dfin) 


dr^ 


with radius R = where tq and a are the radius and diffuseness param¬ 

eters, respectively. (32 is the quadrupole deformation parameter. In Eq. (12), 
we have included only the lowest-order term in the deformation parameter of 
the deformed Woods-Saxon potential ( see, e.g.. Ref. (82|). This is an approx¬ 
imation. However, this should be sufficient for our purpose of illustrating the 
role of projectile deformation effects on the breakup cross sections. 

The radial wavefunction corresponding to the potential 14n(ri) should be 
obtained by solving the coupled equation 


dr\ 


+ ‘^[E - 145(n)] \ Uimin) 


2/r 


i' 


(13) 


Therefore, the radial wave functions obtained from Eq. (13), correspond¬ 
ing to a given I will have an admixture of wave functions corresponding to 
other I values of the same parity. Thus this wave function can be quite 
different from that of the spherical Woods-saxon potential. However, if com¬ 
ponents of the admixed states of higher £ are quite weak, then the pure states 
of lowest I can become dominant. In such a situation, one can use the solu¬ 
tions of the spherical Woods-Saxon potential corresponding to a single £ for 
the wave function uifri) in Eq. (11). Indeed, it has been shown in Ref. M 


that as the binding energies approach zero, the lowest ^ components become 
dominant in the neutron orbits of the realistic deformed potential irrespective 
of the size of the deformation. In this work, we have made the approximation 
of taking U(^{ri) as the state of a given single £ value that is the solution of the 
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Schrodinger equation with spherical Woods-Saxon potential. In any case, if 
the spectroscopic factors of shell model calculations are used for a particular 
state, the wave functions obtained in a spherical basis for that state should 
already include the admixture of different £ states. 

We would like to point out here that only the structure part of the ampli¬ 
tude given by Eq. (11) is affected by the deformation in the interaction Vfen(i'i) 
- the dynamical part of it remains the same as it would be in no-deformation 
case. With the deformation effects introduced through Eq. (12), analytic ex¬ 
pressions can be written for the structure part of the amplitude in Eq. (11). 
Let us dehne 



We can write (see. Ref. 


64j for details of this derivation). 


(14) 


If = 


X 


^W™i(Q) f rldriji^{Qri)ue{ri) 

/imp 

(It I 


(15) 


where R is dehned as 


Ii = 





(2£ -f-1) (2li 1) 


1/2 


dvr 


( h 2 £ 
V 0 0 0 


h 2 £ 

—rrii 0 m 


.( 16 ) 


with \£ — 2\ < li < \£ + 2\ and mi = m. Notice that there would be a limited 
number of h values to be considered, given that £ is the orbital angular 
momentum of the projectile ground state. In the limit of ^2 = 0, the above 
equation would contain the hrst term in the square bracket [involving the 
spherical potential Wsll’i)] that is reduced precisely in the same form as 
that obtained in Ref. for the case where there is no deformation. 

In Eq. (14), the spherical harmonic lj™^*(Q) (where Q = yqn — cxK) can 
be written in terms of product of two spherical harmonics, one depending on 


q„ and the other depending on K, using Moshinsky’s formula [83 


(iQi)'‘i)r‘(Q) 



2li 1 

2L 


\ 1/2 
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X {k-L rm-ML L ML\h 


rMt 


where 


2li + 1 

2L 

ing from 0 to /i 


( 17 ) 


is the binomial coefficient and L = \/2L + 1 with L vary- 


Therefore, the structure part Eq. (lT5ll . can be evaluated 
and would contain the effect of the deformation of the projectile. 

We once again wish to emphasize the analytic nature of our calculation at 
this point. With the structure part given by Eq. (113. the dynamics part in 
Eq. (11) is still given by the Bremsstrahlung integral, which can be solved 
analytically. 

The triple differential cross section for the reaction is related to reduced 
transition amplitude f3im, as 


d^a 


271 


dEhdilhd^n 




(18) 


£m 


where Va is the a—t relative velocity in the entrance channel and p{Ef,, 
is the phase space factor appropriate to the three-body hnal state. 


3. Results and discussions 

The formalism described in Section 2, has been employed to investi¬ 
gate Coulomb breakup of ^^Mg on a Pb target at the beam energy of 244 
MeV/nucleon. In our analysis the calculated one-neutron removal cross sec¬ 


tions are compared with the corresponding data as reported in Ref. [23 


To calculate the Coulomb breakup amplitude [see, Eq. (11)], we require the 
single-particle wave function U£{r) that describes the core-valence neutron 
relative motion in the ground state of the projectile (for a given neutron- 
core conhguration). As discussed in the previous section, we take this wave 
function as that of a state of a single £ value and obtain it by solving the 
Schrodinger equation with a spherical Woods-Saxon potential with radius 
(ro) and diffuseness (a) parameters of 1.24 fm and 0.62 fm, respectively. The 
depth of this well is adjusted to reproduce the valence neutron separation 
energy corresponding to this state. 

Various observables for the reaction have been obtained by integrating 
the triple differential cross sections [see, Eq. (17)] over appropriate angles 
and energies of the unobserved quantity. For example, the total Coulomb 
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Table 1: Depth (V2,s) of the Woods-Saxon potential well as a function of Sn corresponding 
to neutron removal from the 2p3/2, 2si/2, I/ 7/2 orbitals. The values of parameters tq and 
a are taken to be 1.24 fm and 0.62 fm, respectively in all the cases. 


(MeV) 

VI (2p3/2) 
(MeV) 

K, (2*il/2) 
(MeV) 

K, ( 1 / 7 / 2 ) 

(MeV) 

0.01 

43.97 

24.73 

43.36 

0.05 

44.18 

25.28 

43.43 

0.10 

44.42 

25.72 

43.53 

0.15 

44.64 

26.06 

43.62 

0.20 

44.84 

26.36 

43.71 

0.22 

44.92 

26.47 

43.75 

0.25 

45.04 

26.63 

43.81 

0.30 

45.23 

26.88 

43.90 

0.35 

45.42 

27.11 

43.99 

0.40 

45.60 

27.33 

44.09 

0.45 

45.77 

27.53 

44.18 

0.50 

45.94 

27.73 

44.27 

0.55 

46.11 

27.92 

44.37 

0.60 

46.28 

28.11 

44.46 

0.65 

46.44 

28.29 

44.55 

0.70 

46.60 

28.46 

44.64 


one-nucleon removal cross section for a given ij configuration of the valence 
neutron is obtained by integrating the triple differential cross sections over 
angles and energy of fragment b and angles of the valence neutron n. 

The nuclei in island of inversion are expected to have signihcant compo¬ 
nents of 2p — 2h [i'{sd)~‘^{fpY] neutron intruder conhgurations. Indeed, in 


Ref. 23[ , it has been argued that the valence neutron in is most likely 

to have a spin parity (J^) of 3/2“ that corresponds to the 2 ^ 3/2 orbital. In 
this work, we have considered neutron removal from the 2 ^ 3 / 2 , 2 si /2 and I/ 7/2 
orbitals that correspond to ^'^Mg ground state J'" of 3/2“, 1/2’'', and 7/2“, 
respectively. Since the Sn of the valence neutron in the ^^Mg ground state is 
still uncertain, we show in Table 1, values of the depths of the potential well 
as a function of Sn for all of the three orbitals. 

In Fig. 2, we present the results of our calculations for the pure Coulomb 
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Figure 2: Pure Coulomb total one-neutron removal cross section, in the breakup 

reaction of on a Pb target at 244 MeV/nucleon beam energy as a function of one- 

neutron separation energy Sn obtained with configurations ^®Mg(0''') 0 2 ^ 3 / 2 ^^ (solid line), 
^®Mg( 0 +)(g) 2 si/ 2 J^ (dashed line) and ^®Mg(0+) (g) 1 / 7 / 21 ^ (dotted line) for ^"^Mggs using the 
spectroscopic factors (C'^S') 1.0 in each case. The experimental cross section (taken from 
Ref. [i^) is shown by the shaded band. 


(j_i„ in the breakup reaction of on a Pb target at the beam energy 

of 244 MeV/nucleon as a function of Sn corresponding to the one-neutron 
removal from the 2 ^ 3 / 2 , 2si/2 and I/ 7/2 orbitals. For C^S we have used a 
uniform value of one for each configuration. The shaded band in this figure 
shows the corresponding measured cross section taken from Ref. 23| with 
its width representing the experimental uncertainty. We note that calculated 
cross sections obtained with the ^®Mg(0’'') (8) and ^®Mg(0’'') ® 2 si/ 2 Z^ 

configurations (solid and dashed lines, respectively in Fig. 2) overlap with the 
experimental band in the Sn regions of 0.35±0.06 MeV and 0.50±0.07 MeV, 
respectively. Theoretical cross sections for the 2pi/2 case are almost identical 
to those of the 2 ^ 3/2 case. On the other hand, for the ^®Mg(0’'') 0 1 / 7 / 2 ^ 
conhguration there is no overlap between calculated and experimental cross 
sections anywhere, which excludes the assignment of = 7/2“ to 
Therefore, our results are consistent with the assignment of either of the 3/2“ 
and 1/2+ spin parity to the ground state of ^^Mg with one-neutron separation 
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energies in the ranges as stated above. The Sn deduced in our work for either 
of these configurations is within the range of the evaluated value of 0.16±0.( 
MeV as reported in the most recent nuclear mass tabulation 
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Nevertheless, it should be noted that there is a wide variation in the 
C^S values for these states reported in the literature. For ^®Mg(0’'') ® 
and ^®Mg(0^) ® 2 si/ 2 Z^ conhgurations, while the shell model C'^S' values are 
reported to be 0.31 and 0.001, respectively (iil ]. values extracted from an 
analysis of the ^^Mg breakup data [23| are 0.42 ^q;} 2 and O.dOlg;}®, respectively. 
In the latter work the theoretical cross sections have been computed from 
the eikonal model of Ref. for the C target and from the (semiclassical) 
Coulomb breakup model of Ref. for the Pb target. On the other hand, 
for the ^®Mg(0’'') ® 1 / 7 / 2 ^ conhguration the C’^S is not mentioned in these 
references. For a given neutron-core conhguration Sn extracted from the 
Coulomb breakup data is intimately related to the value of C'^S. Therefore, 
it would be interesting to investigate the C’^S dependence of S'„ extracted in 
our study. 



Figure 3: (a) Sn as a function of C'^S for the same reaction as in Fig. 2 with the ^^Mggs 
configuration of ^®Mg(0+) ® 2 p^/ 2 V, (b) same as in (a) for the ^®Mg(0+) 0 ‘ 2 si/ 2 i' config¬ 
uration 

In Figs. 3(a) and 3(b), we show this correlation for the same reaction as in 
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Fig. 2 for the configurations ^®Mg(0’'') 0 2 ^ 3 / 2 ^^ and ^®Mg(0+) 0 2si/2n of the 
respectively. In these calculations, for each C’^S the corresponding 
Sn is deduced from the region of overlap of the calculated cross section and 
the measured data band as shown in Fig. 2. We see that Sn increases steadily 
with increasing C‘^S. Also the uncertainty in the extracted Sn increases with 
C'^S, because at larger C'^S flatter portions of the calculated cross section 
overlap with the data band that encompasses larger parts of the band. It 
may be mentioned here that in our calculations Sn corresponding to the C'^S 
of 0.42 for the conhguration ^®Mg(0’'') 0 2p3/2Z/, is 0.14 ± 0.03 MeV, which 
is lower than the mean value of Sn (0.22 MeV) obtained in Ref. for the 
same C'^S. It is clear from this figure that for a reliable extraction of Sn from 
the Coulomb breakup studies, it is essential to have accurate knowledge of 
the spectroscopic factors for different conhgurations. 
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Figure 4: (a) cr_in as a function of the deformation parameter (52 in the Coulomb breakup 
of on a Pb target at the beam energy of 244 MeV/nucleon with the configuration 

^®Mg(0+) 0 2^3/21^ for ^"^Mggs. The Sn is taken to be 0.35 MeV with C'^S values being 1.0. 
(b) Same as in Fig. (a) for ^®Mg(0+) ® 2si/2!^ configuration with C'^S and Sn of 1.0 and 
0.50 MeV, respectively. In both (a) and (b) the experimental data (shown by the shaded 
region) are taken from Ref. (^ . 

To substantiate the information extracted from the studies of the one- 
neutron removal cross section, it is desirable to investigate other effects and 
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observables to determine the most reliable configuration of To this 

end, in Figs. 4(a) and 4(b) we investigate the effect of projectile deformation 
on a-in for the reaction studied in Fig. 2 . As discussed earlier, the presence 
of neutrons in nearby degenerate j shells (eg., I/ 7/2 and 2 p 3 j 2 ) in ^'^Mg that 
couple strongly to each other by the quadrupole-quadrupole interaction, can 
lead to the quadrupole deformation of this nucleus. 

In Fig. 4(a), we show our results for (T_i„ as a function of (32 for the 
^®Mg(0+) 02 p 3 / 2 Z/ configuration of ^^Mg^s with C‘^S values of 1.0 correspond¬ 
ing to a Sn of 0.35, which is the mean value of the one-neutron separation 
energies extracted from the comparison of the calculated and experimental 
total cross sections in Fig. 2. For (32 = 0, the cr_i„ is the same as that shown 
in Fig. 2 for the same value of Sn- With increasing (32, the cross section 
increases, and the overlap between calculations and the data band ceases for 
(32 > 0.32. Therefore, our calculations do not support a quadrupole defor¬ 
mation parameter in excess of 0.32 for this state. 

In Fig. 4(b) we show the same results for the ^®Mg(0’'') ® 2 si/ 2 ^ con¬ 
figuration with C'^S and Sn values of 1.0 and 0.50 MeV, respectively. The 
contribution of the deformation term to the cross section is substantially 
low for the s-wave configuration, which results in almost constant a-in as a 
function of (32 as seen in Fig. 4(b). We further note that in contrast to the 
results in Fig. 4(a), the overlap between calculated cross sections and the 
data band exists even for values of (32 as high as 1.2. We have checked that 
the situation remains the same for (32 values even beyond 1.2. This result 
indicates that the s-wave configuration does not provide any constraint on 
the deformation parameter (32 in our calculations. On the other hand, the 
Nilsson model study of Ref. [STJ predicts the (32 parameter of the 2 si /2 state 
to be below 0.3. This, however, does not imply that the s-wave configuration 
is negated for the ground state of ^^Mg in our calculations. It simply does 
not constrain the (32 value for this state. In any case, it is not possible to 
obtain more definite constraints on the configuration of from a single 

measurement as available at present. 

The variation of Sn with the deformation parameter (32 is studied in 
Fig. 5 for the same reaction as in Fig. 4(a) for the configuration 

of ^®Mg(0^) ® 2 p 3 / 2 i^ with C'^S value of one. Several authors have argued 
that the deformation can lead to the enhancement of binding energies in the 
island of inversion region nuclei 32|, [2^ [SJ, l85| due to the mixing of 2tuv 


2p — 2h neutron excitations to Ohcu states. We notice in this figures that Sn 
indeed increases with (32- For (32 > 0.70, the Sn value exceeds the upper limit 
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of that evaluated in Ref. j6^. Therefore, for the p-wave conhguration of the 
the deformation parameter remains reasonable even for maximum 
predicted Sn- On the other hand, with the ^®Mg(0’'') ® 2si/2i^ conhguration, 
the Sn remains unchanged with [52, which is obvious from Fig. 4(b). 



Figure 5: Sn deduced from the comparison of our calculations with the experimental data 
as a function of the parameter /32 for the same reaction as in Fig. 4, corresponding to the 
^®Mg(0+) 0 configuration of ^"^Mggs with C'^S = 1.0. 

The investigation of more exclusive observables in the Coulomb breakup 
reactions of the projectile provides signihcant advantages in the understand¬ 
ing of its ground state structure. In Fig. 6, we show the ^®Mg—u relative 
energy spectra in Coulomb breakup of ^^Mg on a Pb target at the beam en¬ 
ergy of 244 MeV/nucleon as a function of the ^®Mg—n relative energy {Erei) 
and [32 simultaneously. The conhguration is ^®Mg(0’'') (8) 2p3/2Z^ with 

C^S and Sn values of 1.0 and 0.35 MeV, respectively. We note that the 
height of the peak depends on the value of (32. The position of the peak 
in this spectrum is dependent on the conhguration of the projectile ground 
state, which is made more explicit in the next hgure. 

In Fig. 7, we show the relative energy spectra {da/dErei) as a function 
of Enel for the same reaction as in Fig. 2 for two diherent conhgurations of 
s^Mg ground state as indicated in the hgure. Since, the peak position of 
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Figure 6: Relative energy spectra for the Coulomb breakup of ^^Mg =3/2 ) on a Pb 
target at 244 MeV/nucleon beam energy, calculated for different values of /32 with 0^5” = 


1.0 and Sn = 0.35 MeV. 


da/dErei is known to depend on the value of Sn , we have used 

the same values of Sn and C^S (0.35 MeV and 1.0, respectively), for the 
two conhgurations. This ensures that differences seen in the relative energy 
differential cross sections of the two conhgurations are attributed solely to 
the differences in the projectile ground state structure. We see that the 
relative energy spectra obtained with two conhgurations show drastically 
diherent behavior as a function of Erei- With the s-wave conhguration, the 
magnitude of the cross section near the peak position is more than 3 times 
larger than that obtained with the p-wave one. Even the peak position of the 
two conhguration are at diherent values of Erei - the p-wave cross sections 
peak at higher Erei as compared to those of the s-wave. 

In view of the results shown in Fig. 7, measurements of the relative energy 
spectra in the breakup reactions of ^^Mg would be of great help in reducing 
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Figure 7: Comparison of relative energy spectra for the Coulomb breakup of ^^Mg on Pb 
target at 244 MeV/nucleon beam energy corresponding to two different possible ground 
state configurations, ^®Mg( 0 +)( 8 ) 2 p 3 / 2 J^ (dashed line) and ^®Mg( 0 +)( 8 ) 2 si/ 2 J^ (solid line). 
The values of S'„ and C'^S' are 0.35 MeV and 1.0, respectively for both the configurations. 
No deformation of the projectile has been included in these calculations. 

the uncertainty in its ground state configuration and also in its one-neutron 
separation energy. Fixing of these quantities will lead to a better under¬ 
standing of the quadrupole deformation of this nucleus, which also affects 
the height of the peak in the relative energy spectra. 

In Fig. 8, we show the parallel momentum distribution (PMD) of the core 
fragment ^®Mg in the Coulomb breakup reaction ^^Mg -|- Pb —>■ ^®Mg -|- n -|- 
Pb at the beam energy of 244 MeV/nucleon. was assumed to have the 

^®Mg(0’'') ® 2p3/2i^ configuration with Sn and C'^S being 0.35 MeV and 1.0, 
respectively. Results are shown for several values of the (32 parameter. We 
note that the magnitude of the cross section near the peak position is quite 
sensitive to the (32 value. Therefore, measurement of this observable is a use¬ 
ful tool for putting constraints on the degree of the quadrupole deformation 
in ^^Mg. 

We note from Table 2 that the full width at half maximum (FWHM) of 
the PMD are almost the same for (32 > 0.30 ( 44.0 MeV/c). Even for (32 = 
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Figure 8: Parallel momentum distribution of fragment in the Coulomb breakup of 

^^Mg on Pb target at 244 MeV/nucleon beam energy for the J’^ = 3/2“ configuration of 
^^Mggs with Sn = 0.35 MeV and C'^S of 1.0. 


0.0, the FWHM is only about 15% larger than its value at higher P 2 - This 
is very close to the FWHM of the PMD of the core fragment seen in the 
breakup reactions of the established low mass halo nuclei like ^^Be and 
Therefore, the p-wave = 3/2“) ground state of ^’^Mg is highly likely to 
have a halo structure. 

In Fig. 9, we show a detailed dependence of FWHM of the PMD on the 
one-neutron separation energy, Sn, for various values of the deformation pa¬ 
rameter /52. The reactions is the same as that studied in Fig. 8 with the 
same value of C'^S. We note that regardless of the value of (32 the FWHM 
increases with increasing Sn- This is expected because with increasing bind¬ 
ing energy the neutron orbits tend to become more and more like those of 
the nuclei away from the drip line. Furthermore, for most values of Sn, the 
(32 dependence of the FWHM of the PMD is similar to that shown in Table 
2 . 

In Fig. 10(a) we show the double differential cross section (Pa/dEndSln 
as a function of the neutron energy for three neutron angles between 1° — 3°. 
The conhguration of ^^Mg^^, and C^S and Sn values were the same as those in 
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Table 2: Full width at half maximum of the parallel momentum distribution of ^®Mg, 
obtained in Coulomb breakup of on a Pb target at the beam energy of 244 

MeV/nucleon. The projectile ground state corresponds to the configuration of ^®Mg(0+)(8) 
2p^/2i’ with Sn and C'^S values of 0.35 MeV and 1.0, respectively. 


(MeV) 

/92 

FWHM (MeV/c) 


0.0 

54.65 


0.1 

50.97 

0.35 

0.2 

48.03 


0.3 

45.82 


0.4 

44.85 


0.5 

44.61 


Fig. 8. No deformation of the projectile was considered in these calculations 
(that is /92 = 0). We see that magnitude of the cross section near the peak 
position reduces with increasing neutron angle. An interesting observation 
is that for all the three angles, the peak occurs near the neutron energy that 
correspond to the beam velocity. This is consistent with the picture that 
fragments move with the beam velocity after the breakup. If the charged 
fragment gets post-accelerated as it leaves the reaction zone (which is the 
case for the breakup reactions of stable nuclei [^, ^), one would expect 
the position of the peak in the neutron spectrum at energies below that 
corresponding to the beam velocity. We do not see this post-acceleration 
effect even if Coulomb effects have been included to all orders in the incoming 
and outgoing channels in our theory. Due to very small binding energies of 
the halo nuclei and the reactions at higher beam energies, the breakup occurs 
at distances much larger than the distance of closest approach, thus the post¬ 
acceleration effects are minimal I 
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The effect of projectile deformation on the cross section d?a/dEndVtn is 
studied in Fig. 10(b) for the same reaction as in Fig. 10(a) for one angle of 
9n = 1°. It is evident from this hgure that magnitude of the cross section 
increases with (32- This is most visible near the peak position. Therefore, 
measurements of the double differential cross sections are expected to provide 
additional information about the deformation of the projectile ground state. 

The angular distributions of neutrons emitted in the projectile breakup 
reactions reflect to a great extent the momentum distribution of the frag- 
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Figure 9: Full width at half maximum of the parallel momentum distribution of ^®Mg, ob¬ 
tained in Coulomb breakup of ^^Mg on a Pb target at the beam energy of 244 MeV/nucleon 
as a function of the one-neutron separation energy Sn and the quadrupole deforma¬ 
tion parameter /32- The projectile ground state corresponds to the configuration of 
^®Mg(0+) 0 2p3/2J^ with a C'^S value of 1.0. 


ments in the ground state of the projectile (see, e.g., Ref. |9l|). Therefore, 


their study is expected to provide further information about the neutron 
halo structure in ^^Mg. In Fig. 11, we show the neutron angular distribu¬ 
tion in the Coulomb breakup reaction ^^Mg + Pb —?• ^®Mg -b n -b Pb at 
the beam energy of 244 MeV/nucleon. The ground state conhguration, the 
Sn and C^S were the same as those in Fig. 10. The results are presented 
for four values of the ^2 parameter. We notice that cross sections drop very 
steeply with increasing neutron angle in the forward directions. The narrow 
angular distributions of neutrons below the grazing angles in the Coulomb 
breakup reactions of ^^Mg reflect the small widths of the parallel momentum 
distribution and hence the large spatial extension of the valence neutron in 
its ground state. The effect of the deformation is signihcant at the forward 
angles [this was already seen in Fig. 10(b)]. 
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Figure 10: Neutron energy-angular distribution for the Coulomb breakup of ^^Mg on a 
Pb target at 244 MeV/nucleon beam energy calculated for Sn = 0.35 MeV and C^S = 
1.0 for the projectile ground state configuration corresponding to J’^ = 3/2“ for (a) (d„) 
at 1° , 2° and 3°, and (b) with different values of P 2 for = 1°, 


4. Summary and Conclusions 

In this paper we have studied the Coulomb breakup reaction -b Pb 
—)■ -b n -b Pb at the beam energy of 244 MeV/nucleon, within the 

framework of the post form hnite range distorted wave Born approximation 
theory that is extended to include the projectile deformation effects. In this 
formalism the transition amplitude is factorized into two parts - one con¬ 
taining the dynamics of the reaction and the another the projectile structure 
informations such as the fragment-fragment interaction and the correspond¬ 
ing wave function in its ground state. Analytic expressions can be written 
for both parts. This formalism opens up a route to perform realistic quan¬ 
tum mechanical calculations for the breakup of neutron-drip line nuclei in 
the medium mass region that can be deformed. 

We calculated the total one-neutron removal cross sections in this 

reaction and compared our results with the corresponding data reported in 
a recent publication in order to determine the conhguration of the ^^Mg 
ground state. The analysis of this single measured cross section already 


25 












Figure 11; Neutron angular distribution for Coulomb breakup of ^^Mg on a Pb target at 
244 MeV/nucleon beam energy. The projectile ground state configuration, Sn and C^S 
were the same as those in Fig. 10. 

rules out the ^®Mg(0+) (g) 1/7/2^^ configuration for the ground state of ^’^Mg. 
However, it does not allow to exclude either of the ^®Mg(0’'') ® ‘2p3i2^ and 
^®Mg(0’'') (8) 2si/2i^ configurations for Assuming a spectroscopic fac¬ 

tor of one, the extracted values of one-neutron separation energies for these 
two configurations are 0.35 ± 0.06 MeV and 0.50 ± 0.07 MeV, respectively. 
However, the deduced S'„ depends on the value of C'^S. Our study shows 
that Sn rises steadily with increasing C'^S. 

In order to gain more insight in the ground state structure of ^^Mg, we 
studied the effect of the projectile deformation on (T_i„. We find that for the 
configuration ^®Mg(0+) 0 2pj,i2V for the ^^Mg ground state, the calculated 
(j_i„ overlaps with the experimental data band for the quadrupole defor¬ 
mation parameter {S2I below 0.32. This is in line with the Nilsson model 
calculations of Ref. (37[ where the (52 for this state is predicted to lie in the 
range 0.30 - 0.34. However, with the ^®Mg(0’'') 0 2si/2Z^ configuration, the 
overlap between calculations and the data occurs for even very large values 
of (52- Thus with this configuration, our calculations are unable to put any 
constraint on deformation parameter (32- 
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We also calculated more exclusive observables such as the core-valence 
neutron relative energy spectra, the energy-angle and the angular distribu¬ 
tions of the emitted neutron and the parallel momentum distribution of the 
core fragment. The position of the peak as well as the magnitude of the 
cross section near the peak of the core-valence neutron relative energy spec¬ 
tra are found to be dependent on the conhguration of the projectile ground 
state as well as on its deformation. Similar trend was also observed in the 
parallel momentum distribution of the core fragment. The FWHM of this 
distribution are found to be of the same order of magnitude as those seen in 
the breakup of established light halo nuclei. This conhrms that ^^Mg ground 
state has a halo structure. The angular distribution of the emitted neutrons 
is strongly forward peaked and the cross sections in the forward directions, 
are dependent on the projectile deformation. Thus, we identihed the observ¬ 
ables that are more critically dependent on the ground state structure of the 
projectile. Therefore, our study is expected to provide motivation for future 
experiments on breakup reactions of the neutron rich medium mass nuclei. 

In calculations of the breakup reactions of nuclei at higher beam energies, 
relativistic effects could play a role 
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Our theory is essentially non- 
relativistic in nature. Nevertheless, we have seen in Ref. that this theory 
is able to reproduce well the data on the excitation energy spectra and the 
total electromagnetic one-neutron removal cross section in the breakup reac¬ 
tion of on a Pb target at even higher beam energy of 422 MeV/nucleon. 
In the present study we did check that inclusion of relativistic effects at the 
kinematics level does not have any signihcant effect on the cross sections. A 
fully relativistic dynamical quantal theory of breakup reactions is still not 
available. 
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